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Abstract 

We develop a novel method in classifying the multipartite entanglement state 
of 2 X X under stochastic local operation and classical communication. In 
this method, all inequivalent classes of true entangled state can be assorted directly 
without knowing the classification information of lower dimension ones for any given 
dimension N. It also gives a nature explanation for the non-local parameters re- 
maining in the entanglement classes while A > 4 . 



1 Introduction 

Entanglement is at the heart of the quantum information theory (QIT) and is now 
thought as a physical resource to realize quantum information tasks, such as quantum 
cryptography [HIS], superdense coding [3lll|, and quantum computation [5], etc. More- 
over, the study of entanglement may also improve our knowledge about quantum non- 
locality [B] . Among it the investigation on the classification of multipartite entanglement 
is of particular interest in QIT. According to QIT, two quantum states can be employed 
to carry on the same task while they are thought to be equivalent in the meaning of 
mutually convertible under Stochastic Local Operations and Classical Communication 
(SLOCC) [7j. 

Nevertheless, in practice the classification of multipartite entanglement in high di- 
mension in the Hilbert space is generally mathematically difficult [8]. It was found that 
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the matrix decomposition method [9l [TO] keeps to be a useful tool as in the two-partite 
case. A widely adopted philosophy in dealing with this issue is first to classify the state 
in lower dimension (or less partite) and then extend to the higher dimension [TT], [12] (or 
more partite [131 cases in an inductive way. However, nontrivial aspect emerges as 
the dimension increases, i.e. some non-local parameters may nest in the entangled states 
[T5t [16]. In recent years, investigations on the classification of 2 x M x states were per- 
formed [Ill[12], where M and N are dimensions of two partites in three-partite entangled 
states. Based on the "range criterion", an iterated method was introduced to determine 
all classes of true entangled states of the 2 x M x N system in Refs. [TTl [12] . In this 
scenario the entanglement classes of high dimensional states can be obtained through the 
low dimensional ones. That is, first generate all the possible entanglement classes under 
invertible local operator (ILO) by the classification information of lower dimensional ones, 
then use the "range criterion" to find out the inequivalent classes of true entanglement 
among all the possible entangled classes, which tends to be a formidable task with the 
increase of dimensions. The main trait of this scenario is that the lower dimensional 
entanglement classes are prerequisite for the follow-up classification. As mentioned in 
Ref. [T2] the classification of entangled state of 2 x M x A^ becomes more subtle when 
M = N. In this case the permutations of the two N-dimensional partites may be assorted 
into different classes. 

In this work, we present a straightforward method in fully classifying the entanglement 
states in 2 X A^ X A^ configuration. The asymmetry of the two A^-dimensional partites 
shows up in one of the classes. We develop a cubic grid form for the quantum state, in 
which the entangled classes that have continuous parameters can be explained naturally. 
This gives an instructive insight on the entanglement classes of 4 or more partites which 
also have non-local parameters [7] [T^. 

The paper is arranged as follows: after the introduction section, we represent the 
entangled state in a general form in section 2, by which the true entangled state of 
2 X N X N can be expressed in a matrix pair. With the definitions given in section 2, the 
true entangled state of 2 x N x N can be fully classified according to the theorems given 
in section 3. In section 4, two examples on how to employ the novel classification method 
are presented. We show, in a typical case of 2 x 5 x 5, that the non-local parameters 
generally may exist in high dimensional or multi-partite entangled state. The last section 
is remained for a brief summary. 
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2 Representation of the entangled states of 2 x x 



An arbitrary two-partite state in dimension of M times can be expressed in the 
following form 

I^MxTv) = 7ii|ll)+7i2|12) + ---7iiv|lA^) + 

72i|21)+ 722|22) + --- + 72^|2Ar) + 



7Mi|Ml) + 7Af2|M2) + ■ ■ ■ ^mn\MN) 



where 7jj e C, are a series of complex numbers. Eq.([T]) can be further expressed in a 
more compact form 
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Here, ^{ij} denotes the MxN complex matrix, which can also be treated as a tensor of 
rank two, and ® is the symbol of direct product. Obviously, the feature of a M x pure 
state is characterized by the rank-two tensor T^ijj. Similarly, the state of 2 x A^ x A^ may 
be expressed in a traced form 



(3) 



where, ipo is a 2-dimensional vector and ipi^2 are N-dimensional vectors, representing the 
constituent states in Hilbert space. The 2x N x N matrix T^ij^ky, which can also be taken 
rank-three tensor, reads 
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Here, Tji ^ ,„} and V s^2,i,m] are in fact tensors of rank two, which are represented by A x A^ 
complex matrices ri^2- The Fi and r2 stand for the upper and lower A^-line blocks of 
matrix Fjij-fc}, respectively. 
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From ([3]) and (jlj) we know that the information of the state 2 x x is involved in 
the matrix pair (Fi , F2). Therefore, in the aim of classification we can specify the typical 
entangled state by a 'matrix vector', that is 



F2 



where, the symbol 



|^2xArxAr) = 

stands for "is represented by" 



(5) 



Generally speaking, two states are said to be SLOCC equivalent if they are connected 
by ILOs [7j. For instance, in the case of bipartite entanglement, suppose the two partites 
are transformed under two invertible local operators P and Q, i.e. ip'-^ = P'^iIji,iIj2 = Qip2, 
then from Eq.([2]) a SLOCC equivalent state to this bipartite entangled state reads as 



(6) 



From above expression, we see that two SLOCC equivalent states are in fact connected 
only by the transformation of the matrix Fjij} in Eq.(l2]) like 



(7) 



Similarly, two SLOCC equivalent 2 x N x N states \'^2xNxn) |^2xAfxAf) are also 
connected by the transformation of the matrix T^ij^k} in Eq.(jl]), i.e.. 
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Here, T is any invertible 2x2 matrix which acts on -i/^o; P and Q are two invertible N x N 
matrices acting on ipi and ip2, respectively. The transformation of the first partite T reads 
as 



(10) 
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For brevity, as in ([5]), equation (IHl) can be formulated as 
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Figure 1: The pictorial description of 2 x 3 x 3 state, where each node corresponds to a base 
vector. Assigning a coefficient to the base vector, we then obtain the corresponding matrix 
element of Fi or . 



where ip s are suppressed. 

To give a pictorial description of the quantum state, we take 2x3x3 case as an 
example, see Figure [H The matrices Fi and I2 are placed in parallel in rear and front of 
the cubic, respectively. Of the cubic grid, each node corresponds to an element in the 
matrix pair (Fi,I^) in Eq.Q . 

In matrix algebra, every ILO which acts on a given matrix can be decomposed as 
a series of products of elementary operations on the matrix, and there exist three such 
elementary operations [17]. Therefore, the matrices T, P and Q in Eq.(l8]), which connect 
the two equivalent wave functions, can be decomposed as such sequence of elementary 
operations. In the pictorial language, here the three types of elementary operation cor- 
respond to three types of manipulation of the cubic grid: type 1, interchange of two 
surfaces; type 2, multiplication of one surface by a nonzero scalar; type 3, addition of 
a scalar multiple of one surface to another surface. Specifically, T is responsible for the 
elementary operations between front and rear; P for upper and lower, Q for left and right 
surfaces, respectively. 

According to the common definition [7], a true 2 x N x N entangled state requires the 
following conditions 

'^(PVo) = 2, r(p^J = r(p^J = , (12) 

to be true, where pi = TTjk{pijk) being the reduced density matrix. Hereafter, we denote 
r to be the rank of matrix. In Quantum Mechanics, to each state there corresponds a 
unique state operator, the density matrix. In the representation of matrix pair the density 
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matrix(elements) can be expressed as 

where = 1,2; = 1,2, ■ ■ ■ , N; k, k' = 1,2, ■ ■ ■ , N. The corresponding reduced 
density matrix is (taking ip2 as an example) 

= J2{TlT,),,,. (14) 

i 

If Det(p^2) 7^ 0, then we know r{p^^ ) = N. 



3 Classification of the tripartite entangled state 2 x 

N X N 



With the above preparation, we can now proceed to classify the 2 x N x N state. 
Generically, the whole space of the state (ri,l2) can be partitioned into numbers of in- 
equivalent sets by different I and n. 



{T,,T,) = {Cn,l} , 
Cn, I = {(fl, 12)1 

(aiFi + piVi] 



(15) 
(16) 



max and rjnin are maximum 



where aj, A G C and + 7^ 0; / G [0, n] and n G [0, N]; r 
and minimum ranks of matrices for all possible values of a, and From the definition 
of ffTBl) . there is no common element in different sets, i.e. Cnj H Cm,k = Cn,i 5m,n5i,k- 
Obviously, every (entangled) state (ri,I^) must lie in one of the subspaces of the set 
{Cn,i} with certain n and /, which in principle can be determined via the transformation 
of equation (fTTj) . since one can always classify a set by certain rules voluntarily. Here, 
the criteria rinax(airi + f^iVi) = n and rmin(a2ri + (^2^2) = I attribute to the group SL(2, 

fmin, Rank(ri), and Rank(r2) are all equal in magnitude, its function can be fulfilled by 
a unit matrix. 
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where Oa and Oa are invertible operators; the "Rank" denotes the rank operation on 
r matrices in upper and lower blocks separately, there will have no invertible matrix 0/ 
exist, which enables 



Oj 



f2 



r2 



(19) 



in case n ^ n or I ^ I, i.e. (Fi, and (Fi, are SLOCC inequivalent. If the operator 
0/ exists, substituting (fT9!) into (fT7|) we may have 



Rank 
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Fi 

f2 



(20) 



and from f|T8l) one knows that n < n and I > I. Similarly, since Oj is an invertible 
operator (matrix), one may also get n < n and / > /, and hence, n = n and / = /. From 
the above arguments, two states, the matrix pairs (Fi, r2) and (fi, f2), connected via 
invertible operator belong to the same subset Cn, i ■ 

Therefore, the entangled classes in set z} with different n and / are SLOCC in- 
equivalent, and the question of performing a complete classification on entangled states 
now turns to how to classify the entangled states in subset Cn i ■ 



3.1 Classification on set Cni with n = N 



From the definition of Cn,i we know that if (Fi, r2) G Cn,i, then there exist an invert- 
ible operator T which enables 
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(21) 



where F{ has the maximum rank N and I^' has minimum rank I. According to matrix 
algebra, in principle one can find invertible operators P, Q and 5* which further transform 
the (F'^, r2) in the following form 
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[22) 



Here, r{J) = rmin(a2Fi + /52F2) with J a matrix in the Jordan canonical form. A typical 
form of J reads 
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in which JnX^i) is a x rii matrix which has the following form 
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(24) 



In all, for every (ri,l2) G Cn,i, there exists an ILO transformation, like 



I?) 



T®P®Q 



['0 



(25) 



Provided r(J) = A^, we know that the rank of the matrix J' = ( J — Aj-E), with Aj being 
any eigenvalue of J, must be less than that of J's. This conclusion contradicts with the 
proviso of J having the minimum rank, since J' and J are correlated through an invertible 
operator, let's say T', 
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(26) 



From above arguments, one observes that the rank of J is less than A^, i.e. / < — 1. 
In the special case of = 2, this observation agrees with the proposition given in Refs. 
[131 [ISj- From Eq. fl2Sl) cjyj = {E,J) is equivalent to Cnj under the joint invertible 
transformations of T, P, and Q, that means the classification on Cn, i can be simply 
performed on cat, i. 



From Eq. (fT4ll one can find that for the quantum state (matrix pair) in cn,i 

:i + iAii 
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(27) 



with /„(x) = Here, j = ipi,'ip2 and n^, Aj are defined in Eq. (!24l) . This tells 

that r{p^J = r{p^^) = N. When / 7^ we readily have r{p^g) = 2. This means the state 
in cn^i is true entangled 2 x N x N state, while / 7^ 0. Otherwise it will not be a true 
entangled 2 x N x N state, which is beyond our consideration. 



Theorem 1 V {E,J) G Cn,i, the set Cn,i is of the classification of Cnj under SLOCC: 

(i) if two states in C^j are SLOCC equivalent, then they can be transformed into the 
same matrix vector {E, J); 

(ii) matrix vector {E, J) is unique in ctv, i up to a trivial transformation, that is if {E, J') 
is SLOCC equivalent with {E,J), then {E,J') = {E,J+ XE) with A being an arbitrary 
complex number. 
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Proof: 

(i) Suppose there exists the transformation 
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according to equation (125!) 
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(ii) Suppose 
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as noted beneath the Eq. fl23p we have / < A^— 1, and it tells that there are no zero elements 
in the pivot of T'. Then, T' can be decomposed as lower and upper triangular(LU) forms 

m 
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where a, /5, 7, A G C and both matrices on the righthand side are nonsingular. Now 
Eq. (!30i) becomes 
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Here, according to the definition of operator Q in Eqs.lJTTl) and ( l22l) . Q' is applied to 
matrix vector from the right hand side. Thus the operator O2 acts on {E, J) in following 

way 



( E ] ^ P ] ( E ' 
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where J = J{Xi 



iJuiiK)- The Eq. fl32|) now gives two independent equations 

E = P'{aE + (3J)Q' , 
J' = XE + -fP'JQ'. 
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(36) 
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The first one can be reformulated into a different form 

P'J{a + (3Xi)Q' = E 
and from eqs.( !35ll and (!37l) we can furtlier get 
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Q.E.D. 

3.2 Classification on set C„ ; with n = N — 1 

Having classified the C^j, now we proceed to the Cjy^i^i case. For every (Fi,!^) G 
Cn-\,Ii we can find an ILO T which implements the following transformation 

(42) 

where r{T[) = N — 1 and r(l^') = /. Then in principle one can find ILOs Pi and Qi, which 
transform the matrix vector (r{,r2') into the form (A, Fg") where A is a x diagonal 
matrix with N-1 nonzero elements of 1 and one zero, and I^" is another N x N matrix in 
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a specific form. To be more explicit, taking = 6 as an example (but the procedure is 
N independent) the above mentioned procedure tells 
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where A, B-^, c and r are submatrices of partitioned by vertical and horizontal lines; 
X means no constraint on the corresponding element. Note that for the case of N-by-N 
matrix, r2 can also be partitioned into a similar form as (jBD, with the lower right block 
to be still the B^ (see Appendix [Al for details). 

In general, there are four different choices for c and r, i.e., 1) c = 0, r = 0; 2) 
c 7^ 0, r = 0; 3) c = 0, r 7^ 0; 4) c 7^ 0, r 7^ 0. In these cases Fg' can be further simplified 
through a series of elementary operations, e.g. denoted Pk and Qk, which enables 
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Here, the superscripts c and r stand for different choices mentioned in above, and the Fg'' 
read 
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In the case of r2"^'', it has aheady been in the form of 
cases we can repartition the submatrices, hke 
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where n in Bn means the dimension of the matrix. This procedure results in the enlarge- 
ment of the blocks Bn, the shrink of blocks A, and the emergence of new types of off 
diagonal blocks c and r. 

By the same procedure, one can further simplify, enlarge Bn and shrink A, the r2'^^ of 
the forms (14711491) . Finally the I^" in (jUj) may arrive at the form of 
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with different kinds of i?„s, correspondingly. Here, J is the Jordan canonical form of A, 
and D is an invertible matrix. Every I^" matrix has its own specific form of Bn block. 
Note that Bn may be recursively enlarged to be the whole matrix of i.e., n = N. In 
all, for every (ri,r2) G Cn-i,i, there exists an ILO transformation, like 
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Here, the r=^Q ^ ^ ^ = Yli^i Q ~ YliQi^ where Pi,Qi stand for those 
operators in Eqs. (144|) . (l45l) . and (!50|) . From Eq. (!5T!) the subset c^-ij, defined as 

C;v-i,/ = {(A,r)|r(r)=/;r= Q ^; (A,r)eC;v-M}, (52) 

is equivalent to C^^ij under the joint invertible transformations of T, P, and Q, which 
means that the classification on C^^i^ i can be simply performed on c^r_i^ 



Similar as ( [27]) we find 

"^i (-[ \\ \ \2\m 

Det(p,) = nE L _ 1 , f!nrr\^) J-2^'(^")-V0, (53) 
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where j = ipi,4'2, and n^, Aj is defined as the Jordan blocks in Eq. (!50|) . and conclude that 
all the states in cn^h are truly entangled in the form of 2 x x A^. 



Theorem 2 V (A, F) G cn^ij, the set c^^ij is of the classification of C^^ij under 
SLOCC: 

(i) suppose two states in Cn-i,i are SLOCC equivalent, they can then be transformed into 
the same matrix vector (A,F); 

(ii) the matrix vector in ctv-i,/ is unique up to a nonzero classification irrelevant param- 
eter, i.e., provided (A, F') is SLOCC equivalent with (A,F), then (A, F') = (A, kF) in the 
sense of J' equalling to J as in theorem 1 while B'^ being exactly the same as B^. 



Proof: 

(i) According to the property of transitivity in SLOCC transformation, this proposition 
should be true. 



(ii) Suppose 
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we first demonstrate that the three ILOs transformation T', P', Q' can always be fulfilled 
through two ILOs transformations P", Q", that is 



T'®P'®Q'[ ^\ = P"®Q"^ ^ 



P"AQ" 



P"kTQ" 



(55) 



13 



According to the definition of cn-i,i we can write (A, F) in tlie form of direct sums 
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is a square matrix, and its dimension equals to that of The transformation T' 



can be decomposed into the following form 
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according to the nature of direct sum. 

For Eq. (!59|) . from the proof of theorem 1 one can find the following Pj, Qj 
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should exist. Here J + \E is taken to be equivalent to the J while r( J + \E) = r{J) = I. 

For Eq. (p7]) . we take into account the two decomposed operations of T' separately. We 
have 
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with cr = ^ . As shown in Appendix [Bl there exist ILOs Pbu and Qsn which satisfy 

Qb„ ={b„)- (63) 
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p , A' + X'Bn 

-T Rr) I 



And, further more we have 
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a 
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There also exist ILOs and Q'^^ which satisfy (see Appendix | 

A' " 



P' 



B„ + AA' 



. ]Q'Bn 



A' 

B„ 



Thus Eq. fl60|) becomes 
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Br, 



where Pb = P'snPsn and Qb = QbuQeu- By taking Pq = Pj ® and Qq 
we have 
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Then, Eg. (15^ can be expressed as 

A 



P'Pr 



A 



7 
Qo'Q' 



A 

r 



p"Ag 

P"kTQ" 
which is just Eq.flS^. 

Eq. (l68|) corresponds to two equations 

P"Ag" = A 

pvrg" = r' 

The equation A = P"AQ" requires P" and Q" taking the following form 



P" 



S Y 
p 



Q"- 
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s-^ 

X q 



respectively. Here, p and q are arbitrary complex numbers. Note that in order to guarantee 
P" and Q" to be non-singular matrices, p and q can not be zero. 

The canonical form of V in cn-i.i gives further constraints on the patterns of P" and 

they can be generally 



Q" . Of the r' and F, in which the takes the form of (j4 
expressed as 
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(71) 



In fact any elements in set cn-i,i takes the same patten flTTj) . From fITOl) . P" and Q" take 
the forms of 

/ an ai2 ais Vi \ 

021 0-22 0-23 0,24 0-25 1/2 

^31 '^32 ^33 0,34 035 1/3 

O4I CI42 CL43 044 045 1/4 

0-51 0^2 Cf53 0,54 O55 ^5 

\0 p J 
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(72) 



where S = {ajjjsxs = {^jjlsxs- From (ITT!) and constraint F' = P"kTQ", P" and Q" read 
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• (73) 



We notice that if P" and Q" are invertible, then the upper-left sub-matrices {ajj}3x3 and 
I&?'il3x3 should also be invertible, due to the fact that any invertible matrix in the form 

invertible. From (l69l) one may infer that in P" kTQ" only the block {ciij}3x3 acts on vector 
V = {vi,i = 1,2,3} in F of fl7ip . Since no invertible operator can transform a nonzero 
vector into a null one, therefore v = and v ^ determine two ILO inequivalent cases 
for F. Thus we see that if F' and F in cn-i,i satisfy Eq. fl69l) . the identity of their B3 
sub-matrices leads to the identity of their B4 sub-matrices. Or in other words, for two 
matrices F and F', if their sub-matrices B^ are the same, while their B4 sub-matrices are 
different, then they should be ILO inequivalent, like (jlTj) and 
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The above analysis for is applicable to other sub-matrices Bn with n > 3, e.g., the 
forms of -B4 and B^ in in Eqs. fl47l) - (l49l) . Generally every nonzero element in Bn will 
transform one column or one row of P" and Q" into a unit vector. In the end, under the 
constraint ( |69l) . if two matrices F and F' have the same Bi, then they must possess the 
same -Bj+i. According to theorem 1 we then have 

that is 



' A ^ 




F' 

V J 





(75) 



(Notice that there exists a special case in which B matrix takes up the whole F and then 
the (J59l) does not exist. The elements in the pivot of T' now can have zeros. In this 
situation, the only modification of the above proofing process is by adding two more ILOs 
Pt and Qt which can reverse the flip of the (A', B) induced by the zero elements in the 
pivot of T', see Appendix [Cl for details) 
Q.E.D. 



3.3 Classification on set Cn,i with n = N — i 

The same procedure can be directly applied to the Ciy_2,i case, and so on. Taken here 
again the = 6 case as an example, following the construction procedure in Eq.(jH]) it 
is easy to obtain the canonical form of (A, F) in cq-2j 
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(76) 



(77) 
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In this case, it is obvious that / can not be smaller than 4. Rescale the matrices (176!) and 
( !77|) according to the partition lines we have 



A = E , (7J 



E 
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(79) 



Here F is just like in Eq.l H^ . Thus the classification of set CN-i,i with i > 1 can be 
reduced to the classification of CAr„i^/ in principle. From (jH]) and (177|) we notice that a 
proposition(inequality) of n, / in Cn,i should exist, i.e. 

2(A^-n)</<n, (80) 

when n < N. Eq. (IHOj) is a constraint on the rank of matrix-pair which represents the true 
entangled state of 2 x x A^. 

Now we have fully classified all the truly entangled classes of 2 x A^ x A^ state. A truly 
entangled state of 2 x A^ x A^ must line in one of the sets CN-i,i (or C„j). According to 
(l80ll we can obtain a restriction on the values of n, 

2A^ . ^ 

— < n <N (81) 



when i > 0, and from the arguments above ( ]26l) we know 

n = N , 0<1<N (82) 

when i = 0. From those two theorems proved in this section we know that the mapping 
of CN-i,i ^ CN-i,i determines all the true entanglement classes in CN-i,i- 



4 Examples 

According to above explanation, the classification of the entangled state 2 x N x N 
may be accomplished by repeatedly taking the above introduced procedures. To be more 
specific and for readers convenience, in the following we completely classify the 2x2x2 
and 2x4x4 pure states by using this novel method as examples. 

For N = 2, i.e., three qubits states, there is only one inequivalent set cn=2,i=i (the 
case ci,i does not exist in the three qubits true entanglement state due to proposition 
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(IHOj) ). There are two inequivalent Jordan forms for 2 x 2 matrices with rank one, and thus 
two inequivalent classes in C2,i which correspond to the GHZ and W states separately [7] 
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For = 4, from (I80|) - fl82|) the inequivalent sets are 



Cn,1 — C4^i , C4^2 5 C4,3 



(83) 
(84) 



(85) 
(86) 



Due to (IHTl) . there is no CN-i,i sets in truly entangled states for = 4 when i >2. In the 
case of C4 / all inequivalent classes have the form of 



J 



where 
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10 
10 
10 
1 



^7) 



Hence, we can distinguish them just by virtue of J's pattern. There are two classes in set 
C41, i.e.. 
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six classes in set C42, the 
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and five classes in set 04,3, the 
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In the case of c^j, every class has the form of 



A 

r 



where 



A 



10 
10 
10 




(91) 



From theorem 2 in Section 13.21 we can simply classify the set C3 2 by the pattern of F 
matrix. And, from the measure in constructing matrices and in the same section, 
we find that there is one class in C3 2 







1 

10 



(92) 



and two classes in 03^3 



10 



1 

10 





10 

1 

10 



(93) 



Altogether, there are 16 genuine entanglement classes in 2 x 4 x 4 states, which agrees 
with what obtained in Ref.[T2]. From (193!) one can easily conclude that the permutation 
of the two 4 dimension partites are sorted into different classes, which was noticed in |12j . 

In above examples we enumerate various distinct classes of the 2x2x2 and 2x4x4 
states, each with a representative state. In fact there are still reducible parameters in 
the representative states, the eigenvalues of the Jordan form, e.g. the As in ( I89p . These 
parameters can be sorted into two categories: one with only redundant parameters, which 
can be eliminated out of the state through ILOs; another possesses non-local parameters 
which can not be eliminated through the ILOs and will keep on staying in the entangled 
state as free parameters. For the first case, we take one typical class in set 05^2 as an 
example. The first three diagrams of Fig. ([2]) exhibits the procedure of how the redundant 
parameters being eliminated through elementary operations. Multiplying the vertical or 
horizontal A-B plane of the cubic form in (II) by a factor of ^ (elementary operation 
type 2), we can get the form of diagram (III). The multiplication of the back plane of 
the cubic by factor of (elementary operation type 2) will transform the parameter 
^^j^ from node B to nodes C, D, E, which is represented by the arrow between B and C. 
Then, the parameters can be easily eliminated by elementary operations in three vertical 
planes containing nodes C, D, and E, respectively. 

For the second case, different from the situations shown in first three diagrams, the 
flow of parameters in performing the elementary operations may form loops, like shown 
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Figure 2: The pictorial procedure of eliminating the parameters in entanglement state. The 
plain nodes represent and the solid dark nodes represent 1 if not further specified. (I) represents 
the initial state in cubic form; (II) is a transformed form from (I) by subtracting the front plane 
from the back one; (III) shows the elimination procedure of the parameter on node A. (IV) 
shows the case with two free parameters. 



in diagram IV of Fig. ([2]) for one typical class in set 05^4 as an example. As long as the 
loop(s) appears, the non-local parameter(s) in the entanglement state remains, and vice 
versa. The number of non-local parameters therefore equals to the number of the loops. 
It is worth to mention that although these parameters are free ones, they satisfy certain 
relations in giving out the equivalent classes, like 



(94) 



for 05^4. The situation may become complicated as the number of parameters increase. 
To get a deeper insight of the behavior of those non-local parameters, there still needs a 
lot of work. 
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5 Summary and Conclusions 



In conclusion, we put forward a novel method in classifying the entangled pure states 
of 2 X X A^. A remarkable feature of our method in different from what existed in the 
literature is that it does not need to classify the lower dimension cases first [HI [12] . We 
find in practice that this method in classifying the 2xNxN tri-partite entanglement state 
is quite straightforward. Since the software for Jordan decomposition is available, this 
new method can be applied to the classification of a given state via automatic computer 
calculation, which is very important as the partite dimension A^ tends to be large. Last, 
but not least, in this work a pictorial configuration of the entanglement states on the grids 
is proposed , which gives an intuitive demonstration for the non-local parameters, and is 
efficient in eliminating redundant parameters. 
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Appendix 



A The Construction of B matrix 
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(95) 



(96) 



A direct observation on Eq.(l96l) tells that 766 must be zero, otherwise one can find 
invertible operators and which enable 
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751a; 752a; 753a; 754a; 755x 
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1/ 



(97) 



(9^ 



Given Aj are the eigenvalues of submatrix {'^ijx}bxb and A 7^ — Aj, we will find that 
r{P,j.PiT2QiQx + ^PxPi^iQiQx) = N > N — 1. This contradicts to requirement that the 
maximum rank of (tuFi + ti2r2) is — 1. 
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Let 766 = 0, Eqs. (p5|) and (!96l) become 
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(100) 



Since we are considering the true entanglement of 2 x x states, neither the last 
column nor the last row of the matrix in Eq. fllOOp can be completely zero. There exist 
ILOs P2, Q2 which satisfy the following equations 
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makes 
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(105) 



Here, must be zero also, otherwise to keep the form of (11041) unchanged the matrix in 
Eq. fll05p can be transformed into 



fi'li il2 iL o\ 
7^'i 7^'2 7^'3 ik 
ik ik ik ik 
ik iL2 ik ik 
000001 

\ 10/ 
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Clearly this will lead to the same contradiction as 766 does in fl96l) and fl98l) . Thus Eg. (11051) 
becomes 
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(107) 



Applying the same procedure to the last row as we have performed to the last column. 
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we have 
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(108) 



(109) 



where -P = Hi -^j Q = YliQi sequences of invertible operators Pi and Qi, respec- 
tively. 



B The Superpositions of A' and 



(110) 



Eq. fl63p can be written in the following matrix equations 

PB„BnQB„ = B„ 

Here, for the sake of clarity, we label the A' with subscript n to indicate its dimension. 

Following, we show inductively that the invertible matrices Pb„ and Qb„ can always 
be constructed. 

First, in case n = 1, then A' = 0, -Bi = 0, the construction of invertible operators Pbi, 
Qb^ Eq. ffnOj) is trivial. 



In the case of n = 2, Eq. flllOp becomes 

(00) ^^'^ 
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Pb2 = E and Qb2 of the form ( ^ ^ \ satisfy the above equations. 
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Suppose for arbitrary n (11 101) is true, we show that Pb„+i,Qb„+i can also be con- 
structed, satisfying 



-Pb„+i(A^+1 + XBn+l)QB„+i 
PB„+iBn+lQB„+r 
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In one example of n = 5, A' and B can be expressed as follows 
/ 1 \ 



A' 



5+1 



where B 



5+1 



1 

1 

1 

1 

yooooooy 

r 

^5 



B. 



5+1 



/o 





1 








o\ 
































1 








1 





























1 


^0 








1 





0/ 



(0,1,0,0,0). 



The operator Pb^+i and Qb^+i can be constructed as follows 
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where Y = XrQs^ + Xh.'^. Because the rank ( q ^ 1 unchanged under the 



invertible transformation, we can always find such invertible operator 
satisfies 

1 X A / rgB. A /Or 



1 X 
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It is then easy to verify Eq. flll2l) using Eqs. (lll7llll8p . Note that the Pb^+i and Qb„+i 
constructed above correspond to the case of Eq. dllSp . for the case of fllMp the procedure 
is similar. 

Along the same line, it can also be found that there exist such invertible operators 
^k'Qkthat 

PLKQk = K 

(120) 

C The flip of and Bn 

Using the Eq. flllOp and Eq. fll20p . we show that there exist the following invertible 
matrices Pt, Qt which flip the and Bn, like 



Provided 



and 



PtBnQt y V A' 



PB„(A)(A; + Ai?„)QB„(A) = a:, 

Pb^B^QeM = Bn 

p'^^S^)KQ'bA^) = K 

P'^^{X){Bn + XK)Q'BA>^) = Bn 



(121) 



(122) 



(123) 



it can be found that 



Pt = PBA^)Pi,^{-j)PBM , Qt = QBAm'BA-\)QBM (124) 



enables 

PtKQt \ -XBn 
PtBnQt ) V K 

which is equivalent to (I12ip up to irrelevant coefficients. 



(125) 
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